Non-equilibrium breakdown of quantum Hall state in graphene 
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In this report we experimentally probe the non-equilibrium breakdown of the quantum Hall state 
in monolayer graphene by injecting a high current density (~lA/m). The measured critical currents 
for dissipationless transport in the vicinity of integer filling factors show a dependence on filling 
factor. The breakdown can be understood in terms of inter Landau level (LL) scattering resulting 
from mixing of wavefunctions of different LLs. To further study the effect of transverse electric field, 
we measured the transverse resistance between the v — 2 to v — 6 plateau transition for different 
bias currents and observed an invariant point. 
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The quantum Hall effect (QHE)i has been studied 
extensively in 2D systems and its equilibrium electron- 
transport properties are understood to a large extent. 
The breakdown of the QHE under non-equilibrium con- 
ditions due to a high current density has been studied 
to understand its microscopic origin 2,3 . There has been 
a considerable debate in the literature regarding the de- 
tails of the mechanism of QHE breakdown. The pro- 
posed mechanisms include electron heating 4 , electron- 
phonon scattering-^, inter and intra Landau level (LL) 
scattering 7 ' 8 , percolation of incompressible regions 9 and 
the existence of compressible regions in the bulk 10 . 

Recently QHE has also been observed in graphene 11 
and studied extensively-^. In this paper we probe the 
breakdown of the QHE in graphene by injecting a high 
current density (~ lA/m); this results in a large local 
electric field in the system. The unique band structure of 
graphene near the Fermi energy (E = ±c?i|fc|, where c « 
10 6 ms _1 is the Fermi velocity) gives rise to a 'relativis- 
ts' QHE. In a magnetic field perpendicular to its plane, 
the energy spectrum of graphene splits into unequally 
spaced LLs and is given by E n = sgn(n) ^/ \2Tic 2 'eB\n\) , 
where n is the LL index. When the Fermi level lies be- 
tween two LLs, the longitudinal resistance (R xx ) van- 
ishes and the transverse resistance (R xy ) gets quantized 
to (4| n |^ 2 )e 2 • Further, the presence of transverse elec- 
tric field mixes the electron and hole wavefunctions and 
modifies the en ergy spectru m— i 14 i 15 , which is given by 
E n = sgn(n)y/(2hc 2 eB\n\)(l - (3 2 ) 3 ^ - hc(3k ± , where 
(3 = E/(cB), E is the electric field orthogonal to and 
k± is the wave vector in the direction perpendicular to 
E and B. 

The motivation for exploring the breakdown of QHE 
in graphene is twofold - first, the QHE in graphene is 
very different from the QHE in a 2DEG system. The 
LL energy spectrum of 2DEG is equispaced unlike that 
in graphene. The energy scale set by the cyclotron 
gap (AE U ) in graphene at B =10 T, is much higher 
(-1300 K) than its value for a 2DEG (-20 K) at the 
same magnetic fielcU^. The mechanism of breakdown 
in graphene could be the inter-LL scattering due to 
wavefunction-mixing or possibly entirely different if the 
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FIG. 1: (color online) a) Line plot of the longitudinal resis- 
tance (R xx ) and transverse resistance (R xy ) for a monolayer 
graphene device at 300 mK and 9 T. The inset shows an op- 
tical microscope image. The scale bar corresponds to 6 \ivn. 
Probes S and D were used to current bias the device. By 
using a lock-in technique, probe pairs Vi — V2 and V\ — V3 
were used to measure R xx and R xy respectively. 50 nA of AC 
current at 181 Hz was used for these measurements, b) Col- 
orscale plot of R xy , to show the plateaus of varying magnitude 
clearly, as a function of magnetic field at 300 mK. 



lengthscale for variation of the local electric field due to 
defects is comparable to the magnetic length. In such 
a situation, ((3 > 1), a "collapse" of the LL is possible 
before a longer lengthscale breakdown of QHE 13 . Sec- 
ond, graphene shows room temperature QHE 16 at high 
magnetic field, therefore understanding the breakdown 
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FIG. 2: (color online) Critical current measurements in the 
vicinity of integer filling factors at 9 T and 300 mK. a, c, e, g 
and i show the colorscale plot of R xx as a function of I^c an( l 
V g near filling factors 6, 2, -2, -6, and -10 respectively. Color 
bars indicate the resistance in units of kQ. The white dotted 
lines on the colorscale plot mark the dissipationless region. 
The line plots in b, d, f, h and j show slices along the current 
axis at the gate voltages shown in the figure (top) and two 
slices (bottom) each for the equilibrium (labeled with solid 
circles) and non-equilibrium (solid line) biasing conditions. 
The position of these slices is marked in adjoining colorscale 
plots (marked a, c, e, g and i ). 



mechanism can also be useful for metrological resistance 
standards 17 . In addition, the presence of back gate in 
our devices allows us to change the Fermi level. We can 
hence probe the QH breakdown away from the integer fill- 
ing factors without changing the energy spectrum. With 
these motivations in mind we have probed the breakdown 
of QHE near the filling factors (sgn(n)(4\n\ + 2) ) n i 18 , 
^ = —10,-6,— 2,2, 6iS. To better understand the role of 
high current densities we have modeled the effect of high 
current using a current-injection model 15 ; this allows us 
to explain the experimentally observed transition from 
neighboring filling factors in terms of an invariant point. 
We also provide evidence to show that these experimen- 
tal observations cannot be explained purely on the basis 
of local electronic heating. 

To fabricate monolayer graphene devices in a Hall bar 
geometry, we have followed the mechanical exfoliation 
technique 11,12 on degenerately-doped silicon substrates 
coated with 300 nm thick Si02- We optically locate the 
flakes of graphene and pattern electrodes onto them using 
electron beam lithography. The electrodes are fabricated 
by depositing 10 nm Cr and 50 nm of Au by thermal 
evaporation. The degenerately-doped silicon substrate 
serves as a back-gate to tune the density of carriers by 
applying a voltage V g . The inset of FigOJi shows an op- 
tical image of a Hall bar device. In FigOJi we plot the 
longitudinal resistance (R xx ) and transverse resistance 
(R xy ) at T = 300 mK and B = 9 T. Filling factors, 
unique to monolayer graphene (y = ±2, v = ±6) are 
clearly seen for both types of carriers in FigQJi. The mo- 
bility of the device shown in the inset is measured to be 
^11000 cm 2 (Vs) _1 for both types of carriers at 300 mK; 
by using the semi-classical relation for mean free path 
(Z)j2£ the measured mobility gives I = 70 nm for carrier 
density 3 x 10 11 cm -2 . In Fig{Tb we plot the evolution 
of R xy as a function of V g and magnetic field B. The 
plateaus in R xy corresponding to v = ± 2, ± 6 and —10 
are clearly seen. The Dirac peak for the device is shifted 
to about 13 V due to unintentional doping, which cor- 
responds to a charge inhomogeneity of 6 x 10 11 cm" 2 —. 
From the two probe resistance measurements at B = 8 
T, we find that contact resistance is smaller than 700 Q 
in the QH regime. 

To probe the breakdown of QHE, we biased the source- 
drain probes of our device with DC current (Ipc) a l° n g 
with a small AC current (50 nA) in the minima of R xx 
corresponding to filling factors v = ±2, ±6 and — 1(3^ at 
fixed magnetic field. The AC current remains fixed and 
lf)Q is then varied as a function of V g in the vicinity of 
integer v. The AC signals between the voltage probes 
V\ and V2 and the voltage probes V\ and V3 were moni- 
tored with a lock-in amplifier to record the values of R xx 
and R xy respectively. Fig|2] shows the evolution of the 
R xx minima as function of Ijjq for different filling fac- 
tors. The line plots show slices of the data in equilibrium 
and non-equilibrium biasing conditions. In order to in- 
terpret the breakdown from the measured experimental 
data we define a critical current (I^nt) as the linearly 
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FIG. 3: (color online) Plot of critical current (/^S) for dif- 
ferent filling factors at T = 300 mK and B — 9 T. The inset 
shows the plot of Hall voltage developed at breakdown (Vho,zz, 
labeled with circles) and the cyclotron gaps (AE U ) (solid tri- 
angles) plotted on right axis as a function of v. 



extrapolated value of if^ at zero dissipation 2 ^. 

We point out the qualitative features of our data - first, 
with the increase of ijfx?, the width of the dissipationless 
region reduces, eventually leading to the breakdown. Sec- 
ond, the critical current is v dependent. Third, on either 
side of the integer filling factor the boundary of dissipa- 
tion evolves asymmetrically. In Fig|2ki, c, e, g and i for 
zv = 6,2,— 2,— 6 and —10 respectively we see a non-linear 
evolution of the boundary of dissipation as a function of 
4c and V g . 

Before discussing details of the data, we address the 
possible concerns about local heating of the sample that 
could occur in these studies. We have done control ex- 
periment to confirm that heating is not responsible for 
the key experimental observations; by comparing change 
of resistance with temperature and due to current. We 
injected a large DC current through the voltage probes 
Vi and V2 of the Hall bar geometry device to intention- 
ally heat it locally while we simultaneously measured two 
probe resistance-gate voltage characteristic of the oppo- 
site two probes (V3 and V4) using a lock-in technique. 
Then, we observed the evolution of two probe resistance- 
gate voltage characteristic of the same pair (V3 and V4) 
with temperature while the injected DC current was set 
to zero. By comparing these two data sets (change of re- 
sistance with temperature and due to current), we could 
estimate that the temperature of the device does not in- 
crease beyond 3 K while injecting currents as high as 
10 juA. We have also seen an overall asymmetry in the 
evolution of R xx with the sign of injected current close to 
integer filling factor. Pure thermal effects cannot explain 
this asymmetry. Additionally, cyclotron gaps in graphene 
are large and thermal effects cannot completely suppress 
the electric field induced effect. Also, superior thermal 
conductivity of graphene^l is likely to suppress any local 
thermal hot-spot formation. 

In order to understand the mechanism of the break- 
down we examine the dependance of I^ r f t on v. Fig|3] 
shows the plot of I^ r f t for various filling factors v in- 



dicating that I decreases with \v\. The inset shows 
the plot of two relevant quantities, Hall voltage, Vuaii — 
I crit x 7^2-? an d AEj, as a function of v. There is a correla- 
tion between Vuaii an d AE^, which can be explained by 
considering inter LL scattering. The origin of the inter- 
LL scattering is likely to be the strong local electric field 
that mixes the electron and hole wavefunction o 13 i 14 i 15 
providing a finite rate for inelastic transitions. The en- 
ergy for the transitions is provided by the transverse elec- 
tric field parallel to the electron trajectory. This results 
in inelastic scattering between LLs leading to a break- 
down of the dissipationless QH state 7 . Similar inter- 
LL scattering mechanisms have been used to explain the 
breakdown of the QHE in a 2DEG, in samples of width 
less than 10 /am and moderate mobility 2 -^^. The elec- 
tric field for inter-LL (Ell) scattering can be estimated 
to be that field where the quasiparticle can pick up an 
energy corresponding to the LL separation within few cy- 
clotron radii (r c ) i.e. cEll ~ AE v /r c ~ 10 6 V/m. This 
is much higher than the experimentally observed electric 
fields. However, Martin et al. 23 found a much shorter 
lengthscale associated with the charge inhomogeneity 
(^150 nm). The presence of a charge inhomogeneity 23 
leads to a strong local electric field and thus can reduce 
the threshold for the breakdown due to inter-LL scatter- 
ing. 

For v = ±2, VhclII matches quite well with AE U= ±2, 
which indicates that the n = LL width is small. How- 
ever for v = ±6, —10, VHaii is smaller than the corre- 
sponding cyclotron gap. This deviation for v = ±6, —10 
can be explained by considering disorder-induced broad- 
ening of n = ±1,-2 LLs. The difference between 
AE V and eVuaii is approximately ~35meV (~415K). 
These observations are consistent with the experiments 
measuring quantum Hall activation gap 24 , which have 
also revealed similar width of these LLs in samples of 
similar mobilities. Additionally, the difference between 
Vho.ii an d AEj, can also be attributed to inhomogeneous 
charge distribution. Considering inhomogeneous charge 
distribution, the critical current is predicted to be fill- 
ing factor and length scale dependent 13 . However, the 
n = level remains protected from the local electric 
field fluctuations^ 3 -*^. The correlation between Vho.ii and 
AE V shows consistency with the breakdown mechanism 
based on this picture too. In addition, our experimental 
finding that there is a non-linear evolution of dissipation 
boundary can possibly be attributed to Hall field induced 
broadening of the extended state band. 

To further explore the effect of transverse electric field, 
due to the high current density, we look at the plateau to 
plateau transition in transverse conductance (cr xy ). Fig HI 
shows the v = 2 to v = 6 plateau transition at T = 300 
mK and B = 10 T for different values of current. As we 
increase the injected current, the transition width starts 
to increase as well. Interestingly, in the transition re- 
gion, all the curves intersect at the filling factor 4 and 
R xx shows a small suppression in peak resistance around 
the same gate voltage (R xx is shown in top- left inset). 
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FIG. 4: (color online) Plot of a xy as a function of V g for 
v = 2 to v — 6 plateau transition at T— 300 mK and B — 
10 T for different values of currents starting from 0.75 jiA 
with an increment of 1.5 fx A. The invariant point at v — 4 
is clearly seen. The top-left inset shows the plot of R xx as a 
function of V g for the same transition at the same values of 
current as indicated in the main plot. The bottom-right inset 
shows the calculated values of a xy as a function of the V g . 
The invariant point at v — 4 is also clearly seen. 

Such an invariant point indicates that as we increase the 
current, the center of the electric field induced broadened 
extended state band does not move with the current. For 
v = 2 to v = 6 plateau transition, the Fermi level crosses 
the four fold degenerate n = 1 LL. It has been shown 
that at very high magnetic field, spin degeneracy can be 
lifted 26 giving rise to an additional plateau at v = 4. We 
speculate the current invariant point at v = 4 and sup- 
pression of R xx at the same time as a precursor of Zee- 



man splitting. To understand our data quantitatively, 
we carried out numerical calculations based on the in- 
jection model of QHE in graphene 15 . This model gives 
transverse conductance from the calculation of local den- 
sity of states. The bottom-right inset in FigJH shows the 
conductance curves calculated for different values of the 
injected current. This model accurately describes the po- 
sition of the current invariant point but fails to explain 
the width of the transition region. One possible reason 
for the failure of this model could be the assumption 
that all states are extended. Further detailed analysis is 
needed to take into account the effect of disorder. 

In summary, we have studied the non-equilibrium 
breakdown of the quantum Hall state in graphene. We 
find that the dissipationless QH state can be suppressed 
due to a high current density, and the corresponding crit- 
ical current decreases with \v\. The correlation between 
VhclII and AE U is consistent with the disorder-induced 
broadening of LLs and inhomogeneous charge distribu- 
tion. The value of Vnaii a t breakdown gives an idea about 
the activation energy. Scanned probe based measure- 
ments on cleaner samples are likely to observe the electric 
field induced "collapse" ((3 > 1) of LLs. We also see a 
current invariant point in the plateau to plateau transi- 
tion and suppression in longitudinal resistance at higher 
current, which can possibly be a sign of spin-degeneracy 
breakdown. 
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